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O : 

Presently, the MAGIC telescope uses a 300 MHz FADC system to sample the transmitted 
and shaped signals from the captured Cherenkov light of air showers. We describe a 
method of Digital Filtering of the FADC samples to extract the charge and the arrival 

■ time of the signal: Since the pulse shape is dominated by the electronic pulse shaper, a 

numerical fit can be applied to the FADC samples taking the noise autocorrelation into 
account. The achievable performance of the digital filter is presented and compared to 

' other signal reconstruction algorithms. 
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in ! 1 Introduction 
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^ ■ The purpose of the MAGIC Telescope pQ is the observation of high energy gamma radiation 
from celestial objects. When the gamma quanta hit the earth atmosphere they initiate a cascade 
of photons, electrons and positrons. The latter radiate short flashes of Cherenkov light which 
can be recorded by a Cherenkov telescope. The FWHM of the pulses is about 2 ns. 

In order to sample this pulse shape with the 300 MSamples/s FADC system |2j, the original 
pulse is folded with a stretching function leading to a FWHM greater than 6 ns. To increase 
^ , the dynamic range of the MAGIC FADCs the signals are split into two branches with gains 
differing by a factor 10. Figure QJi) shows a typical average of identical signals. 

In order to discriminate the small signals from showers in the energy range below 100 GeV 
against the light of the night sky (LONS) the highest possible signal to noise ratio, signal 
reconstruction resolution and a small bias are important. 

Monte Carlo (MC) based simulations predict different time structures for gamma and hadron 
induced shower images as well as for images of single muons [Zj. An accurate arrival time deter- 
mination may therefore improve the separation power of gamma events from the background 
events. Moreover, the timing information may be used in the image cleaning to discriminate 
between pixels whose signal belongs to the shower and pixels which are dominated by randomly 
timed background noise. 



2 Digital Filter 

The goal of the digital filtering method j3JE] is to optimally reconstruct from FADC samples the 
amplitude and arrival time of a signal whose shape is known. Thereby, the noise contributions 
to the amplitude and arrival time reconstruction are minimized. 
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Figure 1: a) Average normalized reconstructed high and low gain pulse shapes and pulse shape imple- 
mented in the MC simulations J^j. b) Noise autocorrelation matrix B for open camera and averaged 
over all pixels (galactic telescope pointing). 

For the digital filtering method, three assumptions have to be made: 

• The normalized signal shape has to be always constant. 

• The noise properties must be constant, especially independent of the signal amplitude. 

• The normalized noise auto-correlation has to be constant. 

Due to the artificial pulse stretching by about 6 ns on the receiver board all three assumptions 
are fullfilled to a good approximation. For a more detailed discussion see jS]. 

Let g(t) be the normalized signal shape, E the signal amplitude and r the time shift between 
the physical signal and the predicted signal shape. Then the time dependence of the signal, 
y(t), is given by y{t) = E ■ g(t — r) + b(t) , where b(t) is the time-dependent noise contribution. 
For small time shifts r (usually smaller than one FADC slice width), the time dependence can 
be linearized. Discrete measurements yi of the signal at times ti (i = l,...,n) have the form 
Hi = E ■ gi — Et ■ c/i + 0{t 2 ) + bi, where g(t) is the time derivative of the signal shape. 

The correlation of the noise contributions at times £j and tj can be expressed in the noise 
autocorrelation matrix B: Bij = (b{bj) — (bi)(bj). Figure |2] shows the noise autocorrelation 
matrix for an open camera. It is dominated by LONS pulses shaped to 6 ns. 

The signal amplitude E, and the product Et of amplitude and time shift, can be esti- 
mated from the given FADC measurements y = (yi, y n ) by minimizing the deviation of the 
measured FADC slice contents from the known pulse shape with respect to the known noise 
auto-correlation: x 2 (E, Et) — (y — Eg — Etcj) t B~ (y — Eg — Etc/) + 0(r 2 ) (in matrix form). 
This leads to the following solution for E and Et: 

-j= rp . n. , , , {g T B' 1 g)B l g-(g T B~ l g)B~ l g 

E = w T amp (t rcl )y + o(T 2 ) with tiffed = ) y T T * :f . . T y > g , (1) 

(g T B g)(g B g) - (g B g) 2 
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Figure 2: a) Amplitude weights wwnp(io) ■ • • Wa.mp{t§) for a window size of 6 FADC slices for the pulse 
shape used in the MC simulations, see text, b) Timing resolution for calibration LED pulses 



-= — t / \ s~\ / 9\ , , . (g T B 1 g)B 1 g - (g T B 1 g)B 1 g 

Et = wl me (t rcl )y + 0(t 2 ) With W timc (trel) = " [ T ?f . ™ g , (2) 

(g T B g)(g B g) - (g B g) 2 

where t re \ is the relative phase between g(t) and the FADC clock. Thus E and Er are 
given by a weighted sum of the discrete measurements yi with the weights for the amplitude, 
^amp(^rei), and time shift, Wtime(^rei), phis 0{j 2 ). To reduce 0{r 2 ) the fit can be iterated using 
g(ti — t — t) and the weights w amp / t ime (t Te \ + r) jU |H]. Figure 121 a) shows the amplitude and 
timing weights for the MC pulse shape. The first weight "Wamp/timcf^o) is plotted as a function 
of t Te i in the range [—0.5, 0.5[ Tadc, the second weight in the range [0.5, 1.5[ Tadc an d so on. 

The expected contributions of the noise to the error of the estimated amplitude and timing 
only depend on the the shape g(t), and the noise auto-correlation B. Analytic expressions can 
be found in references |UIE]- 

3 Performance and Discussion 

Figure Eb) shows the measured timing resolution for different calibration LED pulses as a 
function of the mean reconstructed pulse charge. For signals of 10 photo-electrons the timing 
resolution is as good as 700 ps, for very large signals a timing resolution of about 200 ps can be 
achieved. 

Figure El shows the charge and arrival time resolution as a function of the input pulse height 
for MC simulations (no PMT time spread and no gain fluctuations) assuming an extra-galactic 
background for different signal extraction algorithms. The digital filter yields the best charge 
and timing resolution of the studied algorithms jHJ. 

For known constant signal shapes and noise auto-correlations the digital filter yields the 
best theoretically achievable signal and timing resolution. Due to the pulse shaping of the 
Cherenkov signals the algorithm can be applied to reconstruct their charge and arrival time, 
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Figure 3: Charge and arrival time resolution as a function of the input pulse height for MC simula- 
tions for the Digital Filter, a cubic spline interpolation ( charge= spline integral, time=half maximum 
position) and a sliding window of 6 FADC slices (charge= samples sum, time = pulse barycenter) JSf/. 



although there are some fluctuations of the pulse shape and noise behavior. The digital filter 
reduces the noise contribution to the error of the reconstructed signal. Thus it is possible to 
lower the image cleaning levels and the analysis energy threshold 8J. The timing resolutions is 
as good as a few hundred ps for large signals. 
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